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ABSTRACT 

In this paper , the cooling performance of a refrigerant base nanofluid to cool down electronic chips was 
reported. An alternative base fluid (NH 3 ) is used to substitute the common known base fluid (H 2 0) to enhance the cooling 
capabilities of the nanofluid. Al 2 0 3 -NH 3 with three different volume fractions 1%, 3% and 5%) was examined. To analyze 
the problem , two objective functions namely, the thermal resistance and pumping power, and the Strength Parteo 
Evolutionary Algorithm (SPEA2) was used to optimize the thermal and hydrodynamic performances of the microchannel 
heat sink. The ammonia-base nanofluid (Al 2 0 3 -NH 3 ) outperformed other coolants (SiC-H 2 0, Ti0 2 -H 2 0 , H 2 0 and Al 2 0 3 - 
H 2 0) in pumping power demand (0.144, 0.702, 0.724, 0.94 and 1.015 W for Al 2 0 3 -NH 3 , SiC-H 2 0, Ti0 2 -H 2 0, H 2 0 and 
Al 2 0 3 -H 2 0, respectively) Under the same conditions. Furthermore, lightest microchannel heat sink was needed when the 
ammonia base nanofluid is employed. Finally, the key results from this study open the door for the possibility of 
considering this coolant for excessive heat producer electronic chips. 
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NOMENCLATURE 

Cp Specific heat, J/kg-K 

D h Hydraulic diameter, m 

/ Friction factor 

H Heat sink height, m 

G Volumetric flow rate, m 3 / s 

H c Channel height, m 

h av Heat transfer coefficient, W/m 2 -K 

k Thermal conductivity, W/m-K 

n Spherical particle parameter 

A p Pressure drop, kPa 

Re Reynolds number 

R Thermal resistance, K/ W 
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W Heat sink width, m 

w c Channel width, jam 

Wall (fin) width, jam 

V m f Velocity inside the microchannel, m/s 

P Pumping power 

Greek Symbols 

a Channel aspect ratio 

P Wall width to channel width ratio 

p Density, kg/m 3 

p Viscosity, kg/m-s 

77 Fin efficiency 

(/) Particle volume fraction (%) 

Subscripts 

c Channel 

/ Fluid (coolant) 

hs Heat sink 

nf Nanofluid 

p Particle 

w Wall 

1. INTRODUCTION 

The large scale production of the Very Large Scale Integrated Circuits (VLSICs) in the mid-1970s, has resulted in 
a huge advance in the capabilities of these integrated circuits. On the other hand, they produce a large amount of heat [1]. 
A microchannel heat sink is deemed to be one of the effective heat dissipation devices used for electronic cooling. 
However, its basic capabilities proved to be insufficient to cool down modern electronic chips using common coolants (Air 
and water). Nanofluids were introduced to elevate the cooling capabilities of the microchannel heat sinks. Lee et al. [2] 
showed that the shape and size of the nanoparticles contribute to the enhancement of the overall performance of the 

nanofluids. Xuan and Li [3] studied the effects of using different nanoparticles on the cooling performance of nanofluids. 

They showed that copper particles have a stronger effect than the metal oxide particles. Xue [4] analytical work showed 
that the carbon nanotubes had the best thermal conductivity compared to other nanofluids used in his work. In the past few 
years, several studies have been conducted to explore the genuine capabilities of using nanofluids in electronic cooling 
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systems. Lee and Mudawar [5] in their landmark study investigated experimentally the effect of using nanofluid in both 
single-phase and two-phase flows. They used A1 2 0 3 -H 2 0 as a cooling fluid. Their experiment showed that using nanofluids 
in two-phase flow state leads to catastrophic results because of passages blockage. As for the single-phase, they have 
reported that the improvement in the thermal performance was on the account of high pumping power requirement. 
Mohammed et al. [6] employed a hybrid analysis scheme to investigate the overall performance of a nanofluid 
(A1 2 0 3 -H 2 0) cooled microchannel heat sink under laminar flow condition. The performance of the system has been tested 
for the different volume fraction of the nanofluid. The best performance was noticed when the volume fraction was (5%). 
Meanwhile, Li and Kleinstreuer [7] tested the microchannel heat sink cooled by CuO-H 2 0 nanofluid with Entropy 
generation minimization method. Volume fractions of (1% and 4%) were investigated under laminar flow condition. Their 
key find was that there is a range of Reynolds number where the system can demonstrate its best performance. Ijam and 
Saidur [8] used a thermal resistance model to analyze mini channel heat sink with two different nanofluids ( SiC-H 2 0 and 
Ti0 2 -H 2 0 ). Turbulent flow condition was assumed and the system’s performance was tested for different volume fractions 
(0.8%, 1.6%, 2.4%, 3.2%, and 4%). The optimized result revealed that Ti0 2 -H 2 0 exceeded SiC-H 2 0 in terms of pumping 
power for a volume fraction of 4%. Mohammed et al. [9] in their review article highlighted the drawbacks of using 
nanofluids. Passage clogging due to sedimentation and high pumping power requirement were among the most series 
issues of using nanofluids as coolants. It can be seen from the aforementioned literature that the vast majority of the 
nanofluids used, had water as a base fluid. This paper reports the optimization of a microchannel heat sink using a new 
nanofluid in an endeavor to overcome the drawback of high pumping power requirement. This new nanofluid uses liquid 
ammonia (NH 3 ) as a base fluid with alumina (A1 2 0 3 ) as nanoparticles. The employment of ammonia was to elevate the 
overall performance of the coolant. The selection of (A1 2 0 3 ) was considered because aluminum is well known for its 
resistant to the corrosive effect of ammonia. The optimization scheme developed by Ahmed et al. [10] to analyze the 
overall thermal and hydrodynamic performances of a rectangular microchannel heat sink will be used in this study. 
Modifications have been made to the analysis model by using effector algorithm namely, the Strength Pareto Evolutionary 
Algorithm (SPEA2). 

2. MATHEMATICAL MODELING 

The schematic diagram of the microchannel heat sink considered in this study is shown in Figure 1. 



W c Ww 

Figure 1: Schematic Drawing of the MicroChannel Heat Sink Model 

The proposed nanofluid namely A1 2 0 3 -NH 3 with a volume fraction of 1%, 3%, and 5% is employed. 
The thermophysical properties of A1 2 0 3 -NH 3 , SiC-H 2 0, Ti0 2 -H 2 0, and A1 2 0 3 -H 2 0 are calculated based on the following 
equations [11-14] (n is assumed to be equal to 3). 


www.tivrc.ors 


SCOPUS Indexed Journal 


editor @tjprc. org 




























572 


Ahmed Mohammed Adham 


Pnf ~ Q- + fi)Pf + ^Pp 

(i) 

1 

M " f ~ (l -^) 2 - 5 Mf 

( 2 ) 

(pCp) nf = (l-0)(pCp) f +0(pCp) p 

(3) 

u _k p +(n-l)k f -(n-l)0(k f -k p ) r 
— - . .. - . .. . . kf 

(4) 


nf k p +(n-l)k f +0(k f -k p ) f 

The properties at 20°C of the base fluids (liquid ammonia and water) and the nanoparticles employed in the 
current study are taken from Lienhard [15] and tabulated in Table 1. In this study, the dimensions and the operating 
conditions employed by Tuckerman and Pease [16] are utilized and listed in Table 2. To facilitate the analysis, the 
following assumptions are made. 

• The properties were assumed to be constant. 

• The fluid was single-phase and incompressible. 

• The flow was considered to be a laminar steady flow. 


Table 1: Thermophysical Properties of the Base Fluid and Nanoparticles at 20°C 


Properties 

nh 3 

h 2 o 

A1 2 0 3 

SiC 

Ti0 2 

0 3 

Density, H kg/ m 

609 

994.2 

3970 

3160 

4157 

Specific heat capacity, ^ J/ kg-°K 

4740 

4178 

765 

675 

710 

Thermal conductivity, ^ W/ m-°K 

0.521 

0.625 

40 

490 

8.4 

Viscosity, ^ N. s/ m 

1.38.10' 4 

7.25.10‘ 4 

—- 

—- 

—- 


Table 2: The Dimensions and Operating Conditions used in the Current Study 


Parameters 

Values 

Heat sink width, ^ m 

lxlO ' 2 

Heat sink length, L m 

ixi(r 

TJ 

Channel height, c m 

320xl0' 6 

Substrate thickness, t m 

213xl0" 6 

Volumetric flow rate, ^ m 3 /s 

4.7x1 O ' 6 


The overall thermal and hydrodynamic performances of the considered system are determined through the thermal 
resistance and pressure drop/pumping power. The total thermal resistance is assessed using the thermal resistance equation 
provided in [ 10 ]. 


L 2 1 + /? | 1 1 + yff | t 

Cp nf p nf Re l+« h av l + 2at] k hs 


The required pumping power to circulate the coolant in the microchannel heat sink system is evaluated using the 
following equations given by [ 10 ]. 


tot fh 


(l + a)L V„ 


2 

mnf 


2 // 


■V 


+ (1.79 - 2.23(—-—) + 0.53( 


1+A 


1 V 


1 + /? 


( 6 ) 
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Pp = Ap tot xG 


(7) 


A closer look at Equations (5-7) shows that these performance functions (thermal resistance and pumping power) 

depend on two variables (^and^) therefore, the system will be treated for optimization purposes as a multi-objective 
function with two design variables. 

2.1. Heat Transfer and Friction Factor Correlations 

The heat transfer and friction factor correlations are crucial factors that contribute to the accuracy and credibility 
of the optimization process. Therefore, the following correlations are selected: 


a 


h= 2.253 + 8.1641^—1 


( 8 ) 


ct +1 ) D, 


Equation (8) is given by Kim and Kim [17] is for liquid flowing inside the rectangular microchannel. 

/ = 0.316 Re -0 25 (9) 


Equation (9) offered by Blasius [18] was selected due to the adequate accuracy of the equation when it was 
compared to the experimental results given in [19]. This equation can be used in both laminar and early stages of turbulent 
flows i.e, 3000 < Re <105. 

3. OPTIMIZATION PROCEDURE 


The nature of the problem under hand possess multi-objective functions in nature. The first objective function is 
thermal resistance with pumping power being the other objective function. Two design variables were involved in the 
current analysis namely the channel aspect ratio ( a ) and the wall width to channel width ratio (ft) with their limits listed in 
Table 3. Accordingly, the optimization scheme employed in [19] is recalled to perform the optimization. The SPEA2 is 
used to generate the Pareto optimal front of the objective functions in which a trade-off ability is provided. This SPEA2 
algorithm is shown to possess some preference on NSGA-II in the microchannel heat sink optimization [20]. A code is 
developed along with SPEA2 toolbox in MATLAB (R2015) to produce denser Pareto Optimal front and hence more 
choices will be offered to the designers to best suit their design requirements. Figure 2 shows the complete steps of the 
optimization process. 


Table 3: Design Variables Limits 
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Figure 2: Optimization Procedure of SPEA2 (Adopted from [19]) 

4. RESULTS AND DISCUSSIONS 

In this section, the overall performance of the nanofluid cooled (A1 2 0 3 -NH 3 ) microchannel heat sink is analyzed 
and optimized using three volume fractions (1 %, 3 %, and 5 %). The thermophysical properties of the nanofluid for 
different nanoparticle volume fractions are listed in Table 4. The main focus is on the pumping power required to achieve 
the minimum thermal resistance i.e., the best thermal performance. Firstly, the effect of different nanoparticle volume 
fractions on the overall performance is tested. Then, the performance of the A1 2 0 3 -NH 3 is compared to other coolants. 
Furthermore, the effect of volumetric flow rate on the pumping power and thermal resistance is analyzed. Finally, the 
effect of different heat sink materials on the overall performance of the system is examined. 


Table 4: Nanofluid Thermophysical Properties 


Properties 

ai 2 o 3 -nh 3 

1 % 

3 % 

5 % 

Pnf k 8 /m3 

642.6 

709.8 

111 

// x 10' 4 N. s/m 

' n f 

1.415 

1.489 

1.56 

Cp nf J/kg-°K 

4494.4 

4073 

3724.5 

k nf W/m-°K 

0.536 

0.567 

0.6 
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4.1. Particle Volume Fraction 

The effects of different volume fractions on the pumping power and thermal resistance for three different 
nanofluids are depicted in Figures 3 and 4. It can clearly be seen that with the increase of the volume fraction, the required 
pumping power is increased while the thermal resistance is decreased. For volume fractions of 0.01, 0.03 and 0.05, the 
associated pumping power and thermal resistance are 0.128, 0.138 and 0.144 W, 0.073, 0.07 and 0.068 K7W, respectively. 
The refrigerant base nanofluid (A1 2 0 3 -NH 3 ) outperformed other water base nanofluids (SiC-H 2 0 and Ti0 2 -H 2 0) in terms 
of pumping power by an appreciable factor. Furthermore, it can be inferred that the change in the volume fraction has a 
stronger effect on the pumping power than the thermal resistance. This behavior is expected because the increase in the 
nanoparticle volume fraction, increases the density of the nanofluid which consequently escalates the pressure drop of the 
system which ultimately increases the pumping power (Equations (6-7)). The overall result of this section is compatible 
with the fact that the higher pumping power is always accompanied by the lowest thermal resistance. 

0.072O -... T -T- T-T-I-T-T- 


0.0715 

^ 0.071 

|o.0705 

| 0.07i 

0) 

| 0.0695f 
S 

JS 

H 0.069 
0.0685 

0.068 


$ A1203-NH3 
0-SiC-H2O 
Ti02-H20 




"'O 


0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 

Volume fraction 

Figure 3: Particle Volume Fraction vs. Thermal Resistance 
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Figure 4: Particle Volume Fraction vs. Pumping Power 
4.2. Performance of the Coolant 


In this section, the performance of the A1 2 0 3 -NH 3 is compared to water-base nanofluids namely, A1 2 0 3 -H 2 0, 
SiC-H 2 0, Ti0 2 -H 2 0, and pure water. The thermophysical properties of the aforementioned Nanofluids (for the volume 
fraction of 5 %) and pure water at 20°C are tabulated in Table 5. The optimized results for the same operating conditions 
(Table 2) are listed in Table 6. The performance of the ammonia base nanofluid with (A1 2 0 3 ) nanoparticle was compared 
with the water base (A1 2 0 3 ) nanofluid and pure water and depicted in Figure 5. This was necessary to show the impact of 
usage of the refrigerant (NH 3 ) compared to water (H 2 0) when they used with the same nanoparticle (A1 2 0 3 ). 
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Table 5: Nanofluids and Pure Water Thermophysical Properties 


Properties 

ai 2 o 3 -h 2 o 

SiC-H 2 0 

Ti0 2 -H 2 0 

h 2 o 

Density, ^ kg/m 3 

1146.79 

1100 

1150 

998.2 

Specific heat capacity, Cp J/ kg-°K 

3590.55 

3670 

3550 

4182 

Thermal conductivity, k W/ m-°K 

0.706 

0.723 

0.703 

0.613 

Viscosity, jU N. s/ m 

0.001128 

0.000823 

0.000823 

0.001003 


Table 6: Optimized Results for Different Nanofluids 


Nanofluid 

Optimized Parameters 

R (K/W) 

P(W) 

a 

p 

A1 2 0 3 -NH 3 

0.068 

0.144 

4.987 

0.015 

ai 2 o 3 -h 2 o 

0.055 

1.015 

4.883 

0.029 

SiC-H 2 0 

0.0.070 

0.71 

4.855 

0.026 

Ti0 2 -H 2 0 

0.0695 

0.72 

4.882 

0.025 

H 2 0 

0.057 

0.94 

4.95 

0.041 


1 

1 

| 

u 

0> 

I® 

On 

I 0 

On 

|o 

P H 


Figure 5: Pumping Power vs. Thermal Resistance for Different Coolants 

It can be inferred from Figure 5 that A1 2 0 3 -NH 3 significantly outperformed the other coolants (A1 2 0 3 -H 2 0 and 
H 2 0) in terms of pumping power with slightly higher thermal resistance. This sacrifice in the thermal resistance can be 
justified by the huge gain in the pumping power which is greatly influences the size of the pump that will be used to pump 
the coolant. The comparison was made between these coolants because it shares the same nanoparticle. Furthermore, the 
ammonia base nanofluid also exceeded other water base nanofluids as it can be seen in Table 6 in terms of pumping power. 
The lower pumping power can be explained by the fact that the ammonia-base nanofluid posses lower density compared to 
the water-base nanofluid and pure water. However, the thermal resistance is higher for ammonia-base nanofluid because its 
properties (specific heat and thermal conductivity) that influence the thermal resistance are smaller compared to other 
coolants. Furthermore, the optimized results from Table 6 shows that using the ammonia-base nano fluid provides thinner 
channel walls and hence lighter heat sink i.e., 0.916xl0" 6 , 1.792xl0' 6 and 2.5xl0' 6 m for the A1 2 0 3 -NH 3 , A1 2 0 3 -H 2 0, and 
H 2 0, respectively. 

4.3. Volumetric Flow Rate 

The effects of different volumetric flow rates of the coolant on the pumping power and thermal resistance are 
plotted in Figure 6. All the properties and dimensions are taken from Table 2. The optimized results are listed in Table 7. 
Looking at the far left of Figure 5 and Table 7, it can clearly be seen that the increase in the volumetric flow rate, the 
pumping power increases, and the thermal resistance decreases. This is due to the increase in the volumetric flow rate 
decreases the capacitive thermal resistance (first term in Equation (5)) which eventually decreases the total thermal 
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resistance. Furthermore, it can be inferred from Table 7 that there is an optimum operating point that where the channel 
wall thickness is at its optimum (minimum) value i.e., for G = 8.6xl0'6 m 3 /s, the channel wall thickness is 0.916xl0‘ 6 m. It 
is preferable from the author point of view that the design should be around this optimum point because it has a moderated 
pumping power demand and produced heat sink is lighter. 



Figure 6: Pumping Power vs. Thermal Resistance for Different Flow Rates 


Table 7: Optimized Results for Different Volumetric Flow Rates 


G xlO' 6 (m 3 /s) 

R (K/W) 

P(W) 

a 

p 

7.6 

0.076 

0.112 

4.883 

0.042 

8.6 

0.068 

0.144 

4.987 

0.015 

9.6 

0.065 

0.180 

4.953 

0.023 

10.6 

0.062 

0.231 

4.978 

0.058 

11.6 

0.059 

0.266 

4.997 

0.021 


4.4. Heat Sink Material 

Total pumping power and thermal resistance for different heat sink materials are plotted in Figure 7. Two 
different heat sink materials namely, aluminum and silicon are selected with thermal conductivity values of 238 and 148 
(W/ m-K) for aluminum and silicon, respectively. Resistivity to ammonia's corrosive effect was another selection criterion 
to select the materials. Higher thermal conductivity material offers lower thermal resistance for the same pumping power 
i.e., 0.066 and 0.068 for aluminum and silicon, respectively. This can be explained by the fact that higher thermal 
conductivity lowers the conductive thermal resistance (third term in Equation (5)) and hence lowers the total thermal 
resistance. Therefore, no additional pumping power is required. 



Figure 7: Pumping Power vs. Thermal Resistance for Different Heat Sink Materials 
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5. CONCLUSIONS 

Thermal and hydrodynamic performance of an ammonia-base nanofluid (A1 2 0 3 -NH 3 ) cooled microchannel heat 
sink was analytically studied. One of the very effective algorithms was employed to obtain highly accurate results namely, 
the Improved Pareto Evolutionary Algorithm (SPEA2). The main focus was on the required pumping power to circulate 
the coolant in the microchannel heat sink system. The effects of different parameters on the overall thermal resistance and 
pumping power were analyzed and discussed. It was found that A1 2 0 3 -NH 3 outperformed (lower than) other coolants 
(A1 2 0 3 -H 2 0, SiC-H 2 0, Ti0 2 -H 2 0, and H 2 0) in terms of the required pumping power by ~ 85%. An optimum point was 
found when different volumetric flow rates were used i.e., G = 8.6xl0' 6 m 3 /s, which resulted in a lighter heat sink. In 
addition, higher thermal conductivity material (aluminum) offers lower thermal resistance compared to silicon for the same 
pumping power i.e., P = 0.144 W. Finally, the noticeable results of this research may motivate other researchers to conduct 
more studies using the same nanofluid. 
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